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Abstract 

A Heisenberg model over the square lattice recently introduced by Si and Abrahams to describe 
local-moment magnetism in the new class of Fe-As high-Tc superconductors is analyzed in the 
classical limit and on a small cluster by exact diagonalization. In the case of spin-1 iron atoms, large 
enough Heisenberg exchange interactions between neighboring spin-1/2 moments on different iron 
3(i orbitals that frustrate true magnetic order lead to hidden magnetic order that violates Hund's 
rule. It accounts for the low ordered magnetic moment observed by elastic neutron diffraction in 
an undoped parent compound to Fe-As superconductors. We predict that low-energy spin-wave 
excitations exist at wavenumbers corresponding to either hidden Neel or hidden ferromagnetic 
order. 
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The recent discovery of a new class of high- Tc superconductors that are notably com- 
posed of iron-arsenic layers has reinvigorated the search for new superconductors [l]. Iron is 
usually detrimental to conventional superconductivity because its magnetic moment breaks 
up Cooper pairs [2]. Electronic conduction is confined primarily to the Fe-As layers in the 
new class of high-T^ superconductors, on the other handSj]. The nature of the magnetic 
moments in the iron atoms that make up the new class of materials may then be critical to 
the superconductivity that these systems display. 

As in the copper-oxide high-Tc superconductors, the new Fe-As superconductors are ob- 
tained by doping stoichiometric parent compounds. Elastic neutron diffraction measure- 
ments on the parent compound LaOFeAs reveal the presence of long-range spin-density 
wave (SDW) order at low temperature that is commensurate with the square lattice of Fe 
atoms that make up each layer ^ . The magnetic moment associated with this coUinear type 
of magnetic order is only a fraction of the Bohr magneton, however. Hund's rule is therefore 
violated in the iron 3d orbitals of this new parent compound for high-Tc superconductivity. 

In this Letter, we identify a route to low ordered magnetic moments in frustrated two- 
dimensional magnets composed of local moments of spin one or higher. It is based on linear 
spin-wave analysis and exact diagonalization of a Heisenberg model over a square lattice 
of iron atoms that includes local Hund's rule couplinglsl. We find that Heisenberg spin 
exchange between different 3d orbitals on neighboring iron atoms leads to either hidden 
Neel or hidden ferromagnetic order if the exchange interaction is sufficiently frustrating. 
This may account for the low moment associated with collinear/SDW order that is observed 
in LaOFeAs[4i]. Low-energy spin-wave excitations are a natural consequence of the hidden 
magnetic order, however. We predict that they collapse to the ground-state energy at 
the respective Neel and ferromagnetic wave numbers. Last, we identify a quantum phase 
transition into hidden order from a more familiar frustrated magnetic groundstate that obeys 
Hund's rule. 

Recent transport measurements indicate that parent compounds to iron-based high-Tc 
superconductors are bad insulators (metals) close to a transition into a metallic (insulating) 
state [sl. Further, classical spin- wave frequencies obtained from near- neighbor Heisenberg 
models can be used to fit the measured spin-wave spectra in such parent compounds [Q]. We 
believe, therefore, that a local-moment description of magentism in parent compounds to 
iron-based high-Tc superconductors is valid at low temperature. Following Si and Abrahams, 
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we then consider a spin-1/2 Hamiltonian that contains near- neighbor Heisenberg exchange 
among local iron moments within isolated layers plus Hund's-rule couplingj^: 

^ = ^^oE[Es.(«)l' + EE'^r's.(a)-s,(/5)+ E E^2"'"s.(«)-s,(/3). (i) 
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Above, Si{a) is the spin operator that acts on the spin-1/2 state of orbital a in the iron 
atom at site i. The latter runs over the square lattice of iron atoms that make up an isolated 
layer. The application of Hund's rule is controlled by a negative local Heisenberg exchange 
constant Jq < 0, while nearest neighbor and next-nearest neighbor Heisenberg exchange 
across the links and is controlled by the tensor exchange constants J"'^ and 

J2'^ , respectively. The strength of the crystal field at each iron atom compared to Hund's 
rule determines the number of orbitals per iron atom above. It can be as low as two for 
strong crystal fields, and as high as four for weak crystal fields jsl. We shall now search for 
groundstates of the J0-J1-J2 model above that exhibit low ordered magnetic moments 
that violate Hund's rule. It is useful to first consider the special case where all nearest- 
neighbor and next-nearest-neighbor exchange coupling constants are equal, respectively: 
Jf = Ji and J2'^ = J2- The Hamiltonian then reduces to H = | Jo Z^i Sj ■ Sj + Ji Sj • 

+ J2 Y.{{i.j)) Sj ■ Sj, where Sj = J2a Observe now that Sj + Sj commutes with Sj ■ Sj, 

and hence that the latter commutes with the Hamiltonian. This means that the total spin 
at a given site i is a good quantum number. The groundstate then obeys Hund's rule in the 
classical limit because states with maximum total spin at a given site minimize both the 
Hund's-rule energy (Jq < 0) and the Heisenberg exchange energies in such a case. 

A violation of Hund's rule will therefore require a strong variation in the Heisenberg 
exchange coupling constants among the different iron orbitals. This can be easily seen if we 
confine ourselves to the case of two Sd-wave orbitals per site and choose off-diagonal exchange 
coupling constants that lead to frustration when Hund's rule is obeyed: Jf'" = = 
while Jf'^ = Ji and J^^'^ = J2 if a 7^ /3, with J2 > 0. In the limit of weak Hund's-rule 
coupling, Jq 0, the classical ground state per orbital is a Neel state for Ji < and a 
ferromagnet for Ji > 0. The spins at a given iron atom are equal and opposite across the 
two orbitals, however. (See fig. [U) The moment associated with any type of true magnetic 
order must therefore vanish! It is important to observe that the hidden magnetic order shown 
in fig. [1] is stabilized by the addition of diagonal Heisenberg exchange coupling constants 
that are opposite in sign to the corresponding off-diagonal ones. 



Extremely low ordered moments are therefore possible at weak enough Hund's rule cou- 
pling, Jo < 0, when off-diagonal frustration exists: J^'^ > at a 7^ /3. The hidden order 
that is responsible for it is antiferromagnetic, showing two sublattices (see fig. [T]). Two spin- 
wave quanta per momentum hh are then expected at huJsw above the groundstate energy 7| ■ 
Here, is the natural frequency, which is obtained by linearizing the dynamical equation 
for classical precession by each spin-1/2 moment, Si{a) = Si{a) x dH/dSi{a). In the simple 
case where all diagonal Heisenberg exchange coupling constants are null, it has the form 
cjs^(k) = with 

Q_ = s\Ji\ (2sin^A:»2 + sJ2 (2sin^fc»2 

n=x,y n=+,— 

= 2sJo + s\Ji\ Y {2dnh'^af + SJ2 Y i'^^osh'^af, 

n=x,y ^ n=+,— 

where k' = k or k — {j^ja^-nja) and where cfn = cos or sin, respectively, in the case of 
hidden ferromagnetic order or hidden Neel order per orbital, at off-diagonal Ji > or 
Ji < 0. Above, k'j^ = k'^^k'y^ a denotes the square lattice constant, and s is the electron 
spin. Figure [1] depicts these spectra at maximum off-diagonal frustration = \Ji\/2. The 
spin-wave velocity is then equal to Csw = sa[2{\Ji\+2J2){Jo + 4:9{Ji)Ji + AJ2)]^^'^. It collapses 
to zero at Jo = — 4( Ji -|- J2) for off-diagonal Ji > and at Jq = — 4J2 for off-diagonal Ji < 0, 
which serve as stability bounds for hidden ferromagnetic and Neel order, respectively. 

The above results indicate that large enough off-diagonal frustration in the Jq- Ji- J2 model 
([1]) induces a quantum phase transition into hidden magnetic order that is unfrustrated, but 
that violates Hund's rule. We shall now study how the low-energy spectrum of states for 
the J0-J1-J2 niodel ([T]), with two spin-1/2 moments per site, evolves with the strength of 
the Hund's rule coupling by applying the Lanczos technique numerically on a 4 by 4 square 
lattice with periodic boundary conditions S^]. As usual, we restrict the Hilbert space to 
states with equal numbers of up and down spins. Next, translational invariance is exploited 
to reduce the Hamiltonian to block diagonal form, with each block labeled by a momen- 
tum quantum number. The allowed wave numbers, (fc^^a, /c^a) are then (0,0), (vr, 0), (vr,7r), 
(7r/2,0), (7r/2,7r/2) and (7r,7r/2), plus their symmetric counterparts. The associated trans- 
lational symmetry reduces the dimension of each block to a little under 38,000,000 states. 
Spin-flip symmetry can then be exploited to further block-diagonalize the Hamiltonian at 
such momenta into two blocks that are respectively even and odd under it. The dimension 
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of each of these subspaces is then a httle under 19,000,000 states. Each term in the Hamilto- 
nian ([1]) permutes these Bloch-wave type states, and the permutations are stored in memory. 
Also, the value of the matrix element for Bloch waves that are composed of configurations 
of spin up and spin down that display absolutely no non-trivial translation invariance is 
stored in memory, while it is calculated otherwise. This speeds up the application of the 
Hamiltonian operator tremendously because the vast majority of Bloch waves lie in the first 
category. The application of the Hamiltonian H on a. given state is accelerated further 
by enabling shared- memory parallel computation through OpenMP directives. Last, we 
use the ARPACK subroutine library to apply the Lanczos technique on the block-diagonal 
Hamiltonian operator just described ol. 

Figures [2] and [3] show how the low-energy spectrum of the J0-J1-J2 model ([T]) evolves with 
the strength of the Hund's rule coupling in the case of maximum off-diagonal frustration: 
J"'" = = J^'", while Jf'^ = Ji and J^''^ = | Ji|/2 for a 7^ /3. Respectively, they correspond 
to ferromagnetic and to antiferromagnetic nearest-neighbor Heisenberg exchange, Ji < 
and Ji > 0. Notice first the coincidence at weak Hund's rule coupling, Jq = 0, between 
the previous linear spin- wave approximation about hidden-order shown in fig. [1] with the 
present exact-diagonalization results. It suggests that long-range hidden magnetic order 
indeed exists. Second, notice that the lowest energy spin-1 excitation is not the first but the 
second excited state at strong Hund's rule coupling. This suggests that a nonzero spin gap 
exists at maximum magnetic frustration. We have checked that the low-energy spectrum of 
the corresponding J1-J2 model at spin s = 1 is very similar by setting J^^) = Ju2)- Both a 
spin- wave analysis at large spin s 10|] and series-expansion studies HI] at s = 1 /2 find a spin 
gap at maximum frustration for the J1-J2 model. The spin gap then likely persists at s = 1, 
which argues in favor of a spin gap in the present off-diagonal case. Both sets of spectra are 
then consistent with a transition from a magnetically frustrated state that shows a spin gap, 
but that obeys Hund's rule, to an unfrustrated hidden-order state that violates Hund's rule. 
Figures H] (A) and [5] (A) display level crossings of the lowest-energy spin excitations, which 
are consistent with such a quantum phase transition. It can be shown that a transition into 
hidden magnetic order from true magnetic order of collinear or of Neel type is expected at 
Jq = — 2| Ji| for maximum frustration, J2 = \ Ji\/2, in the classical limit at large spin s. 

Figures H] (B) and [5] (B) also show the evolution of relevant magnetic order parameters 
with Hund's rule coupling. They further confirm the interpretation that a quantum phase 
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transition into hidden order takes place near Jq = — 2| Ji|. The ordered moment is obtained 
here by computing the autocorrelation (0(k)± • 0(— k)±)o of the order parameter 

0(k)± = ^e^''-[S.(l)±S.(2)] (2) 

i 

over the groundstate. Figure H] (B), in particular, displays how the square of the ordered 
moment for true collinear/SDW order (+) decays once the system transits into hidden order 
at off-diagonal Ji < 0. Figure O (B) displays how the same occurs for true Neel order (x) 
at off-diagonal Ji > 0. The former is notably consistent with the low ordered moment that 
is observed by elastic neutron diffraction in an undoped parent compound to the recently 
discovered Fe-As high-T^ superconductors [4]. It must be emphasized, however, that the low- 
energy spectrum of hidden order contains observable spin-wave excitations with energies 
that collapse to the groundstate energy either at the Neel wave number (7r,7r), or at the 
ferromagnetic wave number (0,0). (See figs. [H-El) 

Recent inelastic neutron scattering measurements on an undoped parent compound of 
Fe-As superconductors find a small spin gap at the collinear/SDW wave number (tt, 0) on 
the other handfo]. Figures H] and [6] (A) are consistent with both a reduced moment for 
collinear/SDW order and a small spin gap at (vr, 0) at the transition into hidden Neel order 
for off-diagonal Ji < 0. The undoped parent compounds of Fe-As superconductors could 
then lie at the transition point into hidden magnetic order. 

In conclusion, we have identified a route to low ordered magnetic moments in the undoped 
parent compounds of the recently discovered Fe-As superconductors that is based on the 
weakening of Hund's rule by frustrating Heisenberg exchange interactions between different 
3(i orbitals in neighboring iron atoms. We predict, however, that the spin- wave excitation 
energy vanishes either at the wave number (vr, vr) or at the wave number (0,0) deep inside 
the respective hidden-order phases where Hund's rule is violated. (Cf. ref. [6|].) 

The authors thank Radi Al Jishi and Zlatko Tesanovic for discussions. Exact diagonal- 
ization of the J0-J1-J2 model ([1]) was carried out on the SGI Altix 4700 at the AFRL DoD 
Supercomputer Resource Center. This work was supported in part by the US Air Force 
Office of Scientific Research under grant no. FA9550-06-1-0479 (JPR) and by the National 
Science Foundation under grant no. DMR-0606566 (EHR). 
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HIDDEN NEEL ORDER PER ORBITAL spin-wave excitation energy 




HIDDEN FERROMAGNETIC ORDER PER ORBITAL spin-wave excitation energy 




FIG. 1: The linear spin- wave spectrum for the Hamiltonian ([T]) is displayed in units of |Ji| at 
off-diagonal Ji < and Ji > respectively, at off-diagonal J2 = \ Ji\/2, and with no Hund's rule 
coupling acting on two orbitals per site. Hereafter, we set 7i — > 1. 
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FIG. 2: Shown is the low-energy spectrum for 4 x 4 x 2 spin-1/2 moments that experience off- 
diagonal ferromagnetic and frustrating Heisenberg exchange at weak and at moderately strong 
Hund's rule coupling. The lowest-energy spin-1 state at momentum (vr, vr) is used as the reference 
for the linear spin-wave approximation. 
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FIG. 3: Shown is the low-energy spectrum for 4 x 4 x 2 spin-1/2 moments that experience off- 
diagonal anti-ferromagnetic and frustrating Heisenberg exchange at weak and at moderately strong 
Hund's rule coupling. The lowest-energy spin-1 state at momentum (0,0) is used as the reference 
for the linear spin- wave approximation. 
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FIG. 4: The evolution of the spin gap with Hund's rule coupling at wave numbers that correspond 
to collinear and to Neel order is shown for off-diagonal Ji < and off-diagonal J2 = |</i|/2. Also 
shown is the autocorrelation (O(k)-t- • 0(— k)-|-)o of the order parameter ([2|) for true (+) and for 
hidden (-) magnetic order as a function of Hund's rule coupling. It is normalized to its value in 
the true ferromagnetic state. 
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FIG. 5: The spin gap at wavenumbers corresponding to ferromagnetic and to Neel order are 
displayed as a function of Hund's rule coupling for off-diagonal Ji > and off-diagonal J2 = Ji/2. 
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FIG. 6: Shown is the low energy spectrum for 4 x 4 x 2 spin-1/2 moments that experience maximum 
ofF-diagonal frustration near the transition point into hidden order. 
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